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Tension-Torsion Interaction in Long Bars 
of Thin Symmetric Cross-Sections 


By Denos C. Gazis 

Research Laboratories 
General Motors Corporation 
Warren, Michigan 


Abstract 


An investigation of the interaction of torsion and tension in a thin 
and long rectangular bar is performed by means of a plate theory of 
large rotations. The re-distribution of the tensile stresses under 
the influence of the torsional couple, and the increase of the tor- 
sional stiffness due to the tensile stress field are estimated. The 
same results are obtained by a superposition method. The latter 
method is used for a similar investigation of long bars of a more 
general thin symmetric cross-section. 


Introduction 


The investigation of problems of equilibrium on the basis of the 
linear theory of elasticity is founded on the following assumptions: 
All strains are small compared to unity; the squares and products 
of the angles of rotation are small compared with the strains; and 
the relation between the stresses and the strains is described by 
Hooke’s law.* Any deviation from these conditions introduces a 
non-linearity in some or all of the equations of elasticity. The non- 
linearity has to be taken into account for a more accurate stress 
malysis of various components. Another application of non-linear 
phenomena is in non-linear springs which have a wide use in in- 
dustry. 

The general non-linear theory of elasticity has been formulated 
ty Murnaghan, (1) Biot, (2) Sternberg, (3) and Rivlin, (4) to mention 
afew of the most important contributions. However, very few prob- 
lems have been treated successfully on the basis of the general 
tn-linear equations. This is because, in addition to the intrinsic 
difficulty due to the non-linearity of these equations, a further 

*It may be ascertained that the first two conditions are satisfied 
when the space derivatives of the displacement components are 
much greater than their squares and products. 
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complication arises from the fact that the geometry of the materia] 
bodies may change appreciably when large deformations or rota- 
tions occur. As a result, either the equations of equilibrium or the 
boundary conditions will differ from the corresponding ones of the 
linear theory, depending on whether one employs the coordinates of 
the undeformed or the deformed body. (5) 

A relatively simple class of non-linear problems is derived if 
one permits only large rotations but retains the restrictions of 
small strains and Hooke’s law. Problems of this type may be clas- 
sified as physically linear and geometrically non-linear. (6) Such 
problems may arise in thin plates and shells which involve large 
relative rotations of “plate or shell elements.” Such is, for ex- 
ample, the case of a thin plate strip which is subjected to torsion, 
with or without the simultaneous application of forces and moments 
producing stresses parallel to the plane of the plate. 

Non-linear phenomena in the case of combined torsion and ten- 
sion of beams with thin-walled cross sections have been investi- 
gated by Goodier (7) who discussed the increase in torsional stiff- 
ness of such beams due to the presence of uniform axial tension. 
Another non-linear problem of a similar type, namely, the finite 
twisting and bending of thin rectangular plates, has been treated by 
Reissner (8) on the basis of the von Karman (9) plate theory of large 
rotations, 

In this paper we consider another non-linear phenomenon as- 
sociated with finite twisting of thin and long rectangular bars, 
namely, the distribution of axial stresses or re-distribution of exist- 
ing uniform axial stresses. First, the von Karman plate equations 
are used for the discussion of this phenomenon, It is found that a 
finite twisting of a long elastic bar with a thin rectangular cross- 
section gives rise to axial stresses in every cross-section. If local 
buckling is prevented these axial stresses have a parabolic distri- 
bution across the width with a maximum tensile stress o',,, at the 
ends, a maximum compressive stress o',,/2 in the middle, and zero 
resultant over the cross-section. If an initial uniform axial stress 
is assumed across the width, the twisting causes a re-distribution 
of the axial stresses. The resulting stress-distribution is equivalent 
to a superposition of the uniform stress field and the parabolic dis- 
tribution of axial stresses developed in the absence of initial axial 
stresses. In any case, due to the presence of axial stresses, the 
effective torsional stiffness of the bar is greater than that given by 
the linear theory. The increase of this stiffness over the corre- 
sponding value of the linear theory is found greater than that given 
by Goodier. (7) This is because the re-distribution of the axial 
stresses increases the tensile stresses farthest from the central 
axis of the bar, where they are most effective in resisting a twisting 
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of the bar. The same results are also obtained by superposition of 
two different strain fields, one corresponding to a simple axial ex- 
tension and one corresponding to a twisting of the bar. The associ- 
ated stresses are computed taking into account the non-linearity 
introduced due to large rotations. This superposition method is 
employed for the more general investigation of the tension-torsion 
interaction in long bars of thin symmetric cross-sections, 


Formulation and Solution of the Problem 





Consider a thin elastic plate of length 2 in the x direction, 
width 2a in the y direction and constant thickness 2h, (Fig. 1). The 
von Karman differential equations of equilibrium are” in the ab- 
sence of distributed surface loads, 








a 2 a°F a°w. °F 2?w a°F a°w 
DV w = 2h (353 ay2* ay? ax?” 2 axay dxay! 
1 
‘ a’w .* 3?w d?w - 
V'F = E((— } 


where E is Young’s modulus, w is the transverse deflection of the 
plate, D = 2Eh*/, (1 - v*), v is Poisson’s ratio, and‘ is the square 
of Laplace’s operator in the plane, x, y. The function F is the well- 
known Airy stress function from which the stress resultants in the 
plane are derived according to 





*See, for example, Ref. 10, p. 342. 
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The remaining stress and moment resultants of the plate are given 
by the well-known equations 


p 0° w o°w 
M, = -D(3y3+" oy”) 
a*w o*w 
=, = -D(;,: Y ox? 
3*w 
Mxy = (1-¥) Daay (3) 
Q. =-D—-(v*w), Q, = -p—2 (y*w) 
- 5 ie ’ y ay 
“ aMx ow aw 
Ry, = Q, ay Ni Nxy ax 
2 0 Mxy aw aw 
Ry =Q,- ax ° Nxy ax y ay 


In Eqs. (3), Rx and Ry are the effective shear forces, i.e., the 
stress resultants in a direction normal to the undeformed middle 
plane at sections x = constant and y = constant, respectively. In the 
linear case, R, and Ry become the well-known Kirchhoff shear 
forces, by omission of the last two terms entering their expression 
in Eqs. (3). It might be recalled at this point that the von Karman 
formulation of the problem of finite deformation of plates utilizes 
the coordinates of the undeformed plate. 

The formulation of the problem of combined tension and torsion 
is completed by specifying the pertinent boundary conditions. It is 
required that the resultant force and moment at x = ¢ f be equal to 
the applied external tension, N, and twisting moment, M. The adop- 
tion of integral boundary conditions does not guarantee that the dis- 
tribution of stresses in a physical problem is exactly the same as 
given by the mathematical solution. However, recalling the well- 
known Saint-Venant principle we may accept the mathematical so- 
lution at distances from the two ends of the plate equal or greater 
than the width 2a. Thus, the mathematical solution obtained on the 
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pasis of the integral boundary conditions is the more reasonable, 
the larger the length 2 { is compared to the width 2a of the plate. 

The remainder of the stress resultants at x = + f, and all the 
stress resultants at y=*a are required to vanish. Thus, the 
boundary conditions are, at x= +f 


a 
J N dy =N 
“a 
a a 
f R, ydy + 2 [y Mxy]_. = M (4) 
-@ 
Nyy = M, = 0, 
adat y=ta 
Ny = Nyy = Ry= My =0. (5) 


Inthe second of Eqs. (4) the integral of the left-hand side gives the 

part of the twisting moment which is resisted by membrane action, 

and the remainder isthe part resisted by bending action of the plate. 
In analogy with the linear case we seek a solution with 


w = kxy. (6) 


Substituting (6) in Eqs. (1) we obtain two differential equations in- 
volving F alone, i.e., 


a°F 
8 x dy alg (7) 
*F = Ek’ 
The first of Eqs. (7) yields 
F = f,(x) + f,(y) . (8) 


Furthermore, the symmetry of the system implies a solution sym- 
metric in y and independent of x. Thus one obtains 
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a 2.4 2 
F = (95) Ek’y* + Cy (9) 


where C is an integration constant. 
Having obtained a system of w and F which satisfies the equ- 
tions of equilibrium, we now compute the stress resultants 


2h [E k*y*/2 + 2C] 


M,, = (l-v) Dk 
(10) 
R, = ky N, 
M, = My = Q, = Q,= R, = 0 
The boundary conditions (4) now yield 
N = 4h[Ek’a’/6 + 2Ca] 
(11) 


M = 4a Dk(1-v) + (2kha’/15) (3Ek’a* + 20C) , 


with the remainder of Eqs. (4) identically satisfied. Introducing the 
maximum deflection at the four corners, w,, we have 


k = wo /af (12) 
and hence Eqs. (11) become 


N = (2Eah/3) (6° + 12y) 


(13) 
M= 4(1-v) D6+ (2/15)Eha’6(36* + 20y) 
where 
y = C/E 
(14) 
5 =w/f. 


If the dimensions of the bar and the external loads are given, the 
two dimensionless quantities y and 6 are determined using Eqs. (13). 
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Their values are then entered in Eqs. (10) for the determination of 
the stress resultants N, and Mx, at every point, 

Some observations can be made immediately from the preceding 
equations, regarding the influence of the twist on the distribution of 
stresses and the influence of the tensile stresses on the effective 
torsional stiffness of the bar. Thus subtracting from N, the average 
tension N/2a one obtains 


N, - N/2a = hE6*(y*/a’ - 1/3) . (15) 


The right-hand side of (15) depends only on the parameter 
5 = w,/f and is independent of N, This means that upon twisting of 
the bar with a couple sufficiently large to produce atwist character- 
ized by 6, a parabolic distribution of stresses is developed given by 
the right-hand side of Eq, (15), which is superimposed linearly to 
any existing uniform distribution of tension, N/2a, 

Now, eliminating y from Eqs, (13) one obtains the torsional 
couple M as a function of 6 and N, namely, 


M = 4(1-v) D6 + (1/3) [Nad + (8/15) Ea*hd*] . (16) 


The first term of the expression for M, Eq. (16), is the classical 
rigidity term due to bending action of a plate.* The remaining terms 
correspond to the torsional rigidity due to membrane action. It is 
seen that the torsional moment is resisted partly by bending action 
and partly by membrane action, as pointed out by Goodier. (7) How- 
ever, the term in Eq. (16) corresponding to membrane action is 
greater than that given by Goodier, due to the re-distribution of 
tensile stresses which is considered here. The difference is due to 
the term non-linear in 6, i.e., the second term in the bracket of 
Eq. (16). It is remarked that the re-distribution of tensile stresses 
and the concomitant increase of torsional rigidity is in line with the 
well-known energy theorems of linear elasticity which describe the 
tendency of an elastic medium under the action of a system of forces 
to attain the minimum possible potential energy. Conditions of sta- 
bility have been tacitly assumed to be satisfied. However, in every 
given problem, the stress field described by Eqs. (10) requires a 
closer investigation. The stress resultant N, may become negative 
for small values of y and N, and instability may arise when the 
hegative values of N, exceed in absolute value some physical limits, 





*See, for example, Ref. 10, p. 47. 
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Equation (16) may be written in dimensionless form as 


m=m + Mm, (17) 


where 


M/4D (1 - v) 


B 
ul 


m, =A,6, m,=A,6' 





(18) 
ma =14 (Ase (Ry 
“=? = i) 
and ¢ is the average strain in the axial direction, i.e., 
e = N/4Eah . (19) 


In Fig. 2 are plotted values of m, and m, versus 6, for a wide 
range of A, andA,. The solution of a given problem may be obtained 
by trial and error, using the graphs of Fig. 2. The value of m is 
determined according to Eq. (18), and an appropriate value of 6 if 
found such that the sum of m, and m, corresponding to this 6 is 
equalto m, The range of Fig. 2 may be extended by considering the 
two distinct families of curves as giving the relationship of m, x 10“ 
versus 6 x 10* and m,x 10™ versus 6 x 10", rather than m;(i=1,2) 
versus 6, where k and n are arbitrary exponents. 

Returning to Eq. (16) we observe that it describes the behavior 
ofa non-linear spring due to the presence of the term which is cubic 
in 6. This term is independent of the axial force N and increases 
quite rapidly with 6. (13) The ratio of this term to the classical 
stiffness term is 


(8/45) Ea*h5* _1+v a® ., 
40 -v) Ds a inl (29) 


From Eq. (20) it is seen that as 6 increases the non-linear contri- 
bution to the torsional stiffness increases quite rapidly compared to 
the classical stiffness which is due to plate action only. 
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parameter 5 . 


At this point a look at the stresses involved is in order. It is 
necessary that the stresses do not exceed some physical yield limits 
if one expects to retain the physical linearity of the problem as ex- 
pressed by Hooke’s law. When the loading of the bar is described 
by Eqs. (10) to (14), the maximum shear stress is 


ee 
Tm * i+vp i=) 0 (21) 


and the maximum tensile stress, at y =? a, is 


g,, = (1/3) E6° + N/4ah (22) 
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The increase of this maximum tensile stress over the initial stress 
is 


a, = (1/3) E6* (23) 


m 


and its ratio to the maximum shear stress 


Tm = 1 
T 


+ 

(~) 6 (24) 
m 3 

As is to be expected, for a given 6 the tensile stress o},in- 
creases with respect toT,,, as one considers wider and thinner plate 
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rame 
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tress 


(23) 
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strips. It also increases in proportion with the deflection pa- 
rameter 6. In Fig. 3 the ratio o;,,/7,, is plotted versus the dimen- 
sionless parameter 


Ye, = (AS) (1 


,/G) (25) 
where G is the shear modulus. Two families of curves are obtained, 
using two additional parameters. One family of curves corresponds 
to a/h = constant and is obtained by eliminating 6 between Eqs. (21) 
and (24), in which case one finds 


m=, (F) (26) 


The other family of curves corresponds to o}, = constant, and hence 
it is a family of hyperbolas in an ordinary scale. Figure 3 may be 
used for the direct reading, by interpolation, of one of the three 
quantities, o’,,,7,,, and a/h, when the other two are given. 

A simple numerical example will demonstrate the magnitude of 
the quantities involved. Assume a steel bar with a cross-section of 
1/8” x 4”, E = 30 x 10° psi and v = 0.3, acted upon by an axial force 
N= 2000 lbs., and a torsional couple M = 300 lbs in. Hence, 


a/h 


" 


16, m = 19.97 x 107°, e¢€ = 1.33 x 107* 


A, = 1.022, A, = 88.75 


“ 


Using Figs. 2 and 3 we obtain 


[<a ia.,. €..¢ ie 
m, = 19.4x 10°, o’ = 3,600 psi 


m,= 0.6x 10°, T 


13,400 psi. 


It is interesting to note that the preceding results can be ob- 
lained by a different approach. Assume that the deformation of the 
bar depicted in Fig. 1 is accomplished in two stages, First, we ap- 
ply a uniform axial elongation characterized by a strain ¢, and then 
a twist described by Eq. (6). To the uniform strain corresponds a 
wiform axial stress 
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0% = Eé€,. (27) 


The twist will be, of course, accompanied by the appropriate warp. 
ing of each cross-section and will generate a shear stress distriby. 
tion as described by the torsional theory of Saint Venant. (11) (For 
a narrow rectangular section simpler torsional formulae are also 
available.) (12) In addition to the shear stresses the twist, when as- 
sociated with large rotations, will result in elongation of some longi- 
tudinal “fibers.” Due to the symmetry we are led to assume that the 
centerline is the “twist axis,” which is defined as the unstrained 
axis of the bar during the second stage of simple twisting.* The ad- 
ditional axial strain due to this twisting will be 


e =~ (2) (28) 


where 6 is given by Eq. (14). 

We now set the resultant force and moment at the two ends, equal 
to the applied loads, remembering that the axial stresses contribute 
to the torque in proportion to the slope dw/d@x of the corresponding 
axial fibers. Thus the resultant force and torque are given by 


a 
N= 2h [ o,(y) dy 
-a 


(29) 


= 
iT) 


A 
M, + 2h(5/a) | % (y) y*dy 


where M; is the Saint Venant torque and ot(y) is the total axial 
stress at a distance y from the centerline. It may be ascertained 
that these considerations lead to the same results as given by Eqs. 
(11) to (16). This approach will be employed in the following section 
for the investigation of tension-torsion interaction in long bars of 
various other thin symmetric cross-sections. 


Consider the symmetric cross-section of Fig. 4. The bar is 
again assumed of length 2f, with an applied axial force N and torque 
M at the two ends, xt f. 


= ae fiher 
*It is well known that in the linear theory of torsion any 11be! 


. . ; } . yn- 
may be considered as the twist axis. This in not the case in nm 


linear torsion where twisting is accompanied by axial strain. 











whe 
ory. 


whe 


give 
Eqs 





GAZIS 
(27) 


warp. 
itribu- 
| (For 
e also 
en as- 
longi- 
rat the 
rained 
he ad- 


(28) 


, equal 
ribute 
onding 


(29) 


axial 
tained 
y Eqs. 
ection 
ars of 


par is 
orque 





TENSION-TORSION INTERACTION 13 




















f | 
os &. “6: 
o| eh(y) | 
| | 
z 
| 
| 
, | 
| | 
: | 
Vy 
Fig. 4. Thin symmetric cross-section, 


A uniform axial strain ¢, is applied first, followed by a twist 
about the centerline given by Eq. (6). The centerline is again the 
twist axis of the bar and remains unstrained during twist. The ad- 
ditional axial strain due to this twist will again be given by Eq. (28). 
We now set the resultant force and moment at the two ends equal to 
the applied external loads and obtain 


N= 2E J" [eg+-5 (Z) | niy)dy 
(30) 


= 
“i 


M,+ 2E(6/a) f° [e,+-5 (2) ] h(y)y*ay 


-a 


where M; is the torsional moment given by the linear torsional the- 
ory, corresponding to an angle of twist per unit length equal to 6/a. 

Equations (30) may be used for the determination of €, and 6 
when the geometry of the cross-section and the applied loads are 
given. We give, as an example, the evaluation of the integrals of 
Eqs. (30) for the cross-section of Fig. 5. 


N = Ea [2e,c, + c,) + 67(c, + 3c.) /6] 


M - M; = (1/15) Edéa*[5edc, + 3c,) + 5%c, + 5c,)] . (31) 



























L 


Fig. 5. Example of a thin symmetric cross-section. 








The value of M; may be obtained by using the simplified formula for 
a thin cross-section, (12) Thus we find 


Mt = (4/3) G5 (c, + c,) (c,7 + ¢,”) . (32) 


When c, = c, = h the above formulae coincide with the ones obtained 
for the rectangular cross-section. 
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Summary 


It has been shown theoretically that the twisting of a thin plate 
strip is accompanied by the development of a parabolic distribution 
of axial stresses, with a maximum tensile stress o},, near the edges 
of the strip and a maximum compressive stress o',,/2 along the 
centerline. In the case of an existing tensile stress before twisting 
the parabolic stress field is superimposed to the existing stresses 
as long as the yield limit is not exceeded anywhere in the plate. The 
preceding discussion presupposes that conditions of stability are 
satisfied throughout the material body, and that the yield limit is 
never exceeded during loading at any point. Under this stipulation 
the final stress field is independent of the history of loading. 

The twisting moment is resisted partly by bending action and 
partly by membrane action. The part resisted by bending action is 
the same as given by the classical bending theory of plates, which 
is also equal tothe moment given by the Saint Venant torsion theory. 
The part resisted by membrane action is greater than that previ- 
ously obtained by Goodier, (7) due to the re-distribution of the ten- 
sile stresses. This is because this re-distribution results in a 
higher concentration of tensile stresses near the edges where they 
are more effective in resisting the twisting of the plate. 

Analogous results were obtained for a more general thin cross- 
section in which the thickness is a symmetric function of the dis- 
tance from the central axis of the bar, y. Again, the twisting 
results in a parabolic distribution of axial stresses which is super- 
imposed to the existing axial stresses. The maximum and minimum 
of this parabolic distribution depend onthe variation of the thickness 
across the width. 
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A Modified Newton-Raphson Process for 
Approximation of Multiple Roots of Polynomials 


By Harold Willis Milnes 
Research Laboratories 
General Motors Corporation 
Warren, Michigan 


ABSTRACT 


A modified form of the familiar Newton-Raphson iterative pro- 
cedure is developed especially adapted to the calculation of multiple 
roots. It overcomes many of the difficulties of the standard method 
arising from truncation and roundoff errors and greatly increases 
the convergence rate. A method for determining the multiplicity of 
agiven root is found and an organized computing scheme is pre- 
sented. 


The Newton-Raphson iterative method for approximating the 
roots of a function is a simple, easily applicable device that may be 
applied in the instance of isolated roots with sufficient assurance of 
success to justify its use by computing machinery. However, it is 
well known that if the root to be ascertained has a multiplicity 
greater than unity, then the method fails miserably unless allow- 
ances are made, in an intelligent manner, to correct for inaccura- 
cies and to increase the number of significant digits being used as 
the calculations proceed. With machine calculation where this is 
impossible without taking undesirably elaborate precautions, it is 
found that, as a rule of thumb, roots of multiplicity greater than two 
are usually quite unobtainable and that those of multiplicity two can 
be obtained only at the sacrifice of considerable accuracy. 

The reason for this undesirable behavior of the process is not 
hard to find. If f(x) = 0 is the function whose roots are to be found, 
then the formula for the successive iterates of Newton-Raphson is: 





Xi+, = Xi - (i=1,2,3,-+-) ’ (1) 


where f’(x) is the first derived function of f(x). Now if f(x) has a 
root xX» of multiplicity m > 0, then: 
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f(x) = (x - xX))™ fm(x) , (2.1) 


f’(x) = (x - x9 )™~* [mfp (x) + (x - Xo) f(x)] , (2.2) 


where f,,,(x) denotes a function that does not vanish at xg. As x, 
tends to x,, both f(x) and f’(x) tend to zero, if m 2 2. Consequently, 
the correction term in (1) assumes an indeterminate form. Although 
this indeterminate form admits an analytic evaluation as a limit. 
this cannot be done in actual practice and the machine simply com- 
putes f(x) and f’(x) with as much accuracy as roundoff and truncation 
error allow. The result is that the quotient f{(x)/f’(x) becomes ran- 
dom, aS xj more and more nearly approaches x,, convergence 
ceases and the sequence {x; , , } begins to oscillate in the vicinity of 
the root in a manner which is especially irritating since it is quite 
impossible to estimate how much significance has been obtained be- 
fore this oscillation begins. In the case m = 2, since the factor 
(x - X,) appears only to the first power in the quotient f(x)/f’(x), a 
practical evaluation of the limit may occur, usually rather inaccu- 
rate. 

In this paper, a modified form of the Newton-Raphson process is 
introduced that is especially adapted to the computation of multiple 
roots. Consider equation (2.2) and assume that a close approxima- 
tion of the root has been obtained already, say, through application 
of (1). Thus |x - x,| is taken to be small, although not within the 
desired degree of accuracy and not yet small enough to produce the 
undesirable effects discussed in the preceding paragraph. Then a 
close approximation to f’(x) is given by: 


f’(x) = m(x - x9)" fm (x) . (3) 


Elimination between (2.1) and (3) yields very simply the following 
approximating scheme: 


Xin. = Xi - M ZR 
itl 1 f’(x; 


f (xi) (4) 
Pe, 


It may be noted that the formula is valid even if m is notan integer. 
An interesting relation between (1) and (4) exists, in that the correc- 
tion term in (4) is just the correction term of (1) multiplied by the 
multiplicity of the root. This has the effect of speeding the rate of 
convergence in the manner indicated in Figure 1 for the case m = 2. 
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Formula (4) is exact in case the function f(x) is of the form {(x) 
= c(x - x,)™ where c is a constant; for then: 


m c (x - x,)™ 
Xi+, = Xi - — = Xo > (5) 
mc (x me FF at 





and would therefore yield the root on the first iteration. Its rate of 
convergence is also extremely rapid whenever f,,(x) is slowly vary- 
ing in the region of the multiple root. This is often true in the in- 
stances when Newton-Raphson gives the greatest trouble. Consider, 
for example, the polynomial 
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f(x) = (x - 2)* (x - .01) (x - 4.01) 

(6.1) 
= x* - 8,02x* + 20.1201x* - 16.2404x + 0.1604 


The following table shows that convergence proceeds quite rapidly, 
The advantage obtained may be compared with a similar calculation 
using (1). 


























_ f(x, ) f"(xj) ~ 2f(x; )/f’ (x; ) 

1.500 -0.9345 3.4849 0.536312 

2.036312 -0.00439 -0.29167 -0.03010 (6.2) 
2.006212 0.000153 0.049701 -0.00615682 

2.00005518 t 0.00000003 -- -- 








Formula (4) unfortunately does not relieve the difficulties in- 
volved in the calculation of f(x) due to truncation and roundoff error 
arising from the factor (x - x9)". They are quite as bad as origi- 
nally. We proceed under the assumption that m is now an integer 
and that f(x) is a polynomial. If x; is any number, it is possible by 
the division algorithm to write: 


f(x) = P(x) + (x - Xi) Em(X) , (7.1) 


where gry (x) and r,,,(x) are polynomials and the degree of r,,,(x) is 
less than m. If f(x) is expanded in Taylor series about x; to obtain 


f(x) = {f(xj) + (x - xj) f’(xj) 


m) 
tenet (x = xy)? £1) (xi) } + (x - xi {-—) (7.2) 
™ fa n-m , 
m- t(m+1) (x5) 4... — Hi) f™) (xi)} , 


then g,r,(x) is given by the expression contained in the second set of 
parentheses. Thus, 


1 
Zm(x;) = = f™) (x:) . (8) 


Again if f(x) is expanded about x» we have 
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m £'™)(x ) (x - X,) (m+1) 
f(x) = (x - xq) a + —_ f (x,) 
(9.1) 
n-m 
Keeds (x mn Xo) f'")(x.)} 
n! 
= (x - xq) f,,(x) . (9.2) 
From this follows 
(m) bi (m+ 
f(x; ) = f (Xo) (x; Xo) f - “"(x,) 
m! m+l1! 

(9.3) 

Hees + + - =™ f'")(x,). 

n! 
But by Taylor’s series expansion of {(™) (x) 
f'™)(x,) ’ f'™)(x,) . (x,-x;) fim+1)(x.) 
(9.4) 
én n-rm 
+ (Xo xi) 17) (x) 
n-m! 


From (9.3), (9.4), and (8) it follows that, to within a first approxi- 
mation in (x; - Xo): 


J 1 (m) 
f r(x; ) ” =€ (x;) E1\X;)- (10) 
m! 


Returning to (2.1), by Leibnitz’s theorem 
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1 al 
f(x.) =m! { (xi - xo) * £,,(x;) 
m - k! 
+ ———— (x; - x9)™"*H #7 (x,) (11.1) 
m-k+1! 
1 (k) 
esse 9 =. (x; ” Xo)™ fm (x;)}, 
m! 


from which, to a first approximation in (x; - x, ): 


! _— 
f(x.) = a. (x; ™ Xo) . fn (x;) , 
* (m - k)! (11.2) 
(k=0,1,---(m-1)) . 
By (10): 

f(k) (x,) = a (x; - x,)™* £(™)(x,) (11.3) 

m! ( m-k ( (11 4) 

: (m - k)! ra Gm(Xi) , 


Thus, we obtain the approximating formulae: 


(m - k)! £'*)(x;),1/(m-k) 


— , (k=0,1, --- (m-1)) . (12.1) 
f'™) (x.) 





“a =| = 


We observe that these formulae are valid approximating relations 
for all values of m ranging from unity to the multiplicity of the 
root xo. This is apparent by assuming that any number of factors 
fewer than the full multiplicity of the root have been separated from 
f r(x) in the formula (2.1) and noting that the above arguments still 
apply. Consequently from (10) with m = k, we obtain 


KI gx (xi) 1/(m-¥) 
Xi+¢. = Xi - [— 
m! Emn(X; ) 





(k=0,1, -+-(m-1)) . (12.2) 
























are 3 





NES 


11,1) 


11.2) 


11.3) 


11,4) 


12.1) 











MULTIPLE ROOTS OF POLYNOMIALS 23 


It is to be remarked that these formulae are ambiguous in case 
(m - k) is an even integer and they should not be used in such in- 
stances. 

Of greater interest are the cases m = 1, k= 0 and m equal to 
the full multiplicity, k =(m-1). The former yields the Newton- 
Raphson relation; the latter, the following: 


f(™m~1) (x; ) g m-1(x;) 
Riz, = Xj - .) we (13) 
f™ (x;) Mm g (x ;) 


Formula (13) avoids the difficulties of truncation and roundoff 
error which are inherent in the Newton-Raphson approximation. 
The calculation of g,,-,(x) and g,,(x) is very conveniently accom- 
plished by a scheme of synthetic division analogous to Horner’s. 
The calculations for the equation 


x* - 5x” + 6x’ + 4x - 8 = (x - 2)" (x + 1) 
are arranged below. 


x = 1.6 


~~ 
ss 
~ 
— 
" 
_ 
i 
uo 
Oo 


4 -8 
1.6 -5.44 0.896 7.8336 





A 
tad 
~— 
V 
— 

' 
w 
ae 


0.56 4.896 -0.1664 
1.6 -2.88 -3.712 


| 1.184 


f(1.6) 





i) 
ws 
~ 
~— 
“ 
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_ 
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g, (1.6) 
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1 -5 6 4 -8 
2.0762 -6.0704 -0.1462 8.0013 

1 -2,9238 -0.0704 3.8538 0.0013 
2.0762 -1.7598 -3.7999 

1 -0.8476 - 1.8302 | 0.0539 
2.0762 2.5508 





1 1.2286 | 0.7200 


2.0762 
1 | 3.3048 
= 2.0762 - ——/2— = 2,0035784 
atlas -3x3.3048 “ 
1 -5 6 4 -8 


2.0035784 -6.0035656 -0.0071440 8.0000000 





1 -2.9964216 -0.0035656 3.9928560 | 0.0000000 
2.0035784 -1.9892392 -3.9927407 





1 -0.9928432 -1.9928048 | 0.0001153 
2.0035784 2.0250872 





1 1.0107352 | 0.0322824 
2.0035784 


1 | 3.0143136 


0.0322824 


3x 3.0143196 ~ 270000085 


x, = 2.0035784 - 


We may remark that in the last stage of the above calculation, 
with a machine of seven digit capacity, an accuracy greater than that 
of x, would be unobtainable and oscillation would undoubtedly occur 
past this point. As with any iterative technique, the final results 
are sometimes dependent upon the first approximation. This would 
occur, for instance, if x, = 1.25 for then x, happens to be a root of 
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g, (x) and (13) would be meaningless. A difficulty of a second type 
would occur if x, were a root of g,(x); the approximating relation 
would again be meaningless. Both situations are easily avoided, 
simply by changing the first approximation, Since the approxima- 
tion is, of necessity, still far from the root, this can be accomplished 
by reverting to the usual Newton-Raphson scheme or one or another 
of the relations (12.2). 

In the above, it has been assumed that the multiplicity m of the 
root x, is known. If this is not the case an analysis of the informa- 
tion obtained during the calculation will generally reveal the exist- 
ence of such a root and its multiplicity. It will be assumed that the 
calculations are begun using the conventional Newton-Raphson tech- 
nique until a reasonable approximation is made which, nevertheless, 
does still not greatly affect the accuracy of the computation of f(x). 
From (1) and (3): 


f(x) _ (x 


- X) 


veal 
a ame (14.1) 


Taking the difference of this expression for two successive approxi- 
mations, we obtain 


f (x;) f (xj,,) Xj - Xi+, 
— ._ "=. — (14.2) 
f'(x;) f'(x;,,) m 





From (8) 


m = (14.3) 





f(x) — £(%i+,) | 


Xi - Xi4,|+ 
. , | Ra Gi (X i+.) 


where absolute values are introduced in case f(x) is a complex func- 
tion. In the instance of the numerical example of the previous para- 
graph, at the first iteration: 


| | - (0.1664) - (0.0013) 
m =| 1.6 - 2.0762 | + ~_——__ —_—_—__—_- 
1,184 0.0539 


(14.4) 
2.8931 = 3 
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If m has been ascertained at the it iteration to have some value 
m;, it may be wise to correct this at the (i + 1)* iteration. Since 
f(xj4,) may now have lost significance, it is advisable to make use 
of (11.2) with k = (mj; - 1) and k = mj to obtain 


f™in») (x) _ (x - x, ) 





, (14.5) 
f™i) (x) (m - m,; + 1) 


from which by (8) 


m = (m, - 1) + m; |x; - Xj,, | 


14,6 
Em;-1(%i) — Emj-1 %iss) “ag 


+_ ieee om see 
Em, (xi) Em, (Xi+i) 





In the previous numerical example with i= 2, mj =3 we obtain 
m = 3.0516982. 


The author wishes to thank Dr. D. C. Gazis whose valued sug- 
gestions have been incorporated in the demonstrations leading to 
formula (12.1). He also wishes to thank Mr.C. R. Lewis for his 
assistance in programming calculations for the I.B.M. 704 computer 
at the Research Laboratories. 
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Optimum Allocation of Time and Equipment 
for Several Fatigue Tests 


By Leonard G. Johnson 

Research Laboratories 

General Motors Technical Center 
Warren, Michigan 


INTRODUCTION 


In a previous paper’ the author has considered the problem of 
minimizing the median cost of a single fatigue test which would 
fail r specimens with a known Weibull life distribution, a manufac- 
turing cost C, per specimen, and a waiting cost Cw per hour. The 
solution to this problem consisted of an optimum number (n) of 
specimens to be run simultaneously (with replacements when fail- 
ures occurred) until r specimens had been failed. 

We are confronted with a more general problem when several 
such tests (which may differ from each other in the applied load, 
material used, environment, etc.) are to be started simultaneously, 
and only a limited number (N) of test stations are available. The 
purpose of this report is to outline a method of determining the 
optimum numbers (n;) of specimens to be run simultaneously for 
each of m fatigue tests whose Weibull life distributions, as well as 
waiting costs per hour and manufacturing costs per specimen, are 
mown. It is further assumed that the number (rj; ) of specimens to 
be failed in each test has been decided upon from considerations of 
the coefficient of variation desired in the mean life of each test. 
The final solution must be such that the total number of specimens 


m 
(2 nj) put on test initially does not exceed the total number (N) of 


available test stations. 


NOTATION 


number of fatigue tests to be started simultaneously 
N = total number of available specimen testing stations 
nj = number of specimens to be run simultaneously in the ie 


fatigue test (with replacements as failures occur) 
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r; = number of specimens to be failed in the i‘ fatigue test 


C.;, = manufacturing cost per specimen in the fo 


fatigue test 
Cw, = waiting cost per hour in the i” fatigue test 
L_s5.= median life in hours forthe specimens of the ith fatigue test 


b; = Weibull slope of the life distribution function of the i” fg. 


tigue test 
. r(l ¥ 
> — 
_ mj + Ind-1,°1 b; (See (19) on page 6 of 
¢ = (-—; ) rj! - 
rj 1n2 1 reference 1) 
T(r; + Bb 
1 
. .»d§ 
) Ls _ probability that a specimen of the - fatigue 
P, (x) = 2 i 2 ; 
' test will survive x hours. 


i 1 


, b; > So 
Zi = (n, 7 a4 ¢ 1) C.. + fj " >; Li5,Cw, mi + 


= median total cost of the it” fatigue test (See (24) on page Tof 
reference 1) 


MINIMIZING THE SUM OF MEDIAN COSTS 
FOR AN UNLIMITED NUMBER OF TESTING STATIONS 


The sum of the median costs of all m fatigue tests is 


Ma 


» (nj +r - ICs +t 


Equating zero to the partial derivatives of this expression with re- 
spect to nj, we obtain 


Solvin 


which 
probl 
also | 


n; to 


and 
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1 1 
— -—-1 
GZ 1 bj bj ; 
an; - Cs; r b; rj 9; Lis; Cw; ny = 0 (2) 











“Ll na cee «ee (3) 


which represents the general solution to the optimal allocation 
problem for an unlimited number of testing stations. Equation (3) 
also represents the general solution to our problem as long as 


= n;2N (4) 


MINIMIZING THE SUM OF MEDIAN COSTS 
WHEN (4) IS NOT SATISFIED BY THE n, OF (3) 


In case (4) is not satisfied by the n; of (3), we must reduce each 
nj; to k;, such that 


> Ki - N , (5) 
in which case the sum of the median test costs becomes 


1 


m b; 
(ki +r - 1) Co. + 5 mH * Oj Ls, Cy, kj (6) 


1=1 


~m 
" 
ums 


In order to apply the method of LaGrangian Multipliers? to (5) 
and (6) we form the expression 


wW=S+a( zd kj - N) 
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or 


w= = (ky +H - 1) Cs, 
1=1 
(7) 
=| 
m - D5 m 
+ E> mig Lis,Cw, ki +A(Z ki - N) 


1=1 1=1 


This expression for y is to be minimized. It contains (m + 1) inde- 
pendent variables, i.e., k,, k2,..., Km, and A. We equate zero to 
the partial derivatives of y with respect to each of these (m + 1) 


variables. Thus, in general, 
+ 1 
ay 1 b; ~— 
ok " Cs; ~ D: rj . $j L 9; Cw, kK; +a 0 (8) 
1 1 
aw m 
> * =x ki - N=0 (9) 





(10) 


Since A = 0 makes kj (0) = nj [see (3)], and since we are assuming 
mm 
= n,; > N (since (4) is not satisfied), it follows that we must finda 
i=1 
> 0 to substitute into (10) in order to meet the condition (5). 
If we start with a sufficiently small ¢ (say 10™), we can try 
m 


A =A,=€ in (10), and if this makes > k; (A,) >N, we try A =A, = 


1=) 
in (10), etc., etc., until we find, after q such steps, that A = Ag = 4 
m 
will make > 


m 1=1 
< N, in which case we take the values 


k; (Aq ) > N, but that A =A = (q+1) € will make 


qri 





' 





as t 
succ 
tron 





nde- 
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+ 1) 


(8) 


(9) 


(10) 
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Kj Oger) — oo (a7 7 ica) (11) 
. J 


as the optimal solution to our allocation problem. This method of 
successive trials can be done automatically on a high speed elec- 
tronic digital computer. 


CONCLUSIONS 


From formulas (3) and (10) we conclude that 


1. In case there is an unlimited number of testing stations, the 
optimum number of specimens to run simultaneously for any 
fatigue test of the collection of m tests is given by formula 
(29) of reference 1. 

(Note: Formula (29) of reference 1 is the same as formula 
(3) in this report.) 


2. Formula (29) of reference 1 may be used for the optimum 
number of specimens tobe run simultaneously in each fatigue 


m 
test provided the n, so obtained are such that > n; <= N. 
1=1 
3. In case the n obtained by the formula (29) of reference 1 are 


m 
such that > nj >N, we must reduce each n; to a smaller 


1=1 
value, kj , by choosing a positive A in (10) which is just large 
m 


enough to make > kj; = N. It can be seen that this is equiva- 
1=1 


lent to adding a constant amount A to each original specimen 
cost (since the quantity C,, is replaced by (C,; + A) in each 
test). 
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A Computer Program for Assessing the Reliability 
of an Assembly Consisting of Components with 
Known Weibull Survivorship Functions 


By Leonard G. Johnson 

Research Laboratories 

General Motors Technical Center 
Warren, Michigan 


INTRODUCTION 


It has been emphasized again ‘and again in reliability studies 
that the probability that an assembly of components in series 
will survive a specified service period is the product of the in- 
dependent probabilities of survival of the distinct components 
making up the assembly. Using this general concept, it is pos- 
sible to derive the complete distribution function for the life of 
an assembly once we know the complete distribution functions 
for the lives of the independent components. In this paper it is 
assumed that the components possess Weibull survivorship 
functions such that the probability of component i surviving x units 
of life (these may be hours, stress cycles, or any other time-like 
units) is pj(x) = EXP [ - (x/0;) Pi] where 0; and B; are known posi- 
live parameters for component i. The best fitting Weibull survivor- 
ship function for an assembly of n components is then determined 
empirically by fitting a resultant function EXP [ - (x/0)® | to the 


n 
product Il p;(x) = EXP | - z (x/0,)®i|, that is, the best simple 
1=1 1=1 
power function (5)? is determined for the sum (5 Pi 4 (5-)®: 
1 


> 


+ 


x 
oo + (5 ) Bn in an interval of x values sufficiently large to include 
n 


any Service periods which would usually be of concernto the designer. 


General Method of Attack 


The method by which the resultant best fitting Weibull survivor- 
ship function for an assembly of n components is determined is to 
first obtain the resultant survivorship function for components (1) 
and (2) only. We denote this resultant by (1,2). Next, we determine 
the resultant of the pair (1,2) and (3), which is equivalent to the re- 
sultant of the trio (1), (2), (3). We denote this resultant by (1,2,3). 
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Then the resultant of the pair (1,2,3) and (4) is determined, which js 
the same as the resultant of the quartet (1), (2), (3), (4). Continuing 
in this manner we finally end up with the resultant of the pair 
(1,2,3, ..., (m-1)) and (n), that is, the resultant of the n-tuple (1), 
(2), (3), .... (n), which is (1,2,3,4, ...., n). 

The great advantage of building up a resultant survivorship 
function in this fashion is the fact that we need only design a mathe. 
matical formulation for the resultant of a simple pair of compo- 
nents, and having done this, simply allow the computer to repeat the 
same operations on the latest resultant and the next component of 
the assembly, until every component has been taken into consider. 
ation. The process of determining the resultant of two components 
shall hereafter be known as TWO-COMPONENT SYNTHESIS. 


Two-Component Synthesis 





We are given that pcomponent (1) has a survivorship function 

1) (x) = EXP[( - x/®, y5 1] and that component (2) has a survivorship 
ted P2 (x) = EXP [( - "a 52]. The resultant survivorship 
function is then py, 2) (x) = p(x) p(x) = EXP [ - (& @, )B: 
~ (x/8, )B2}. We want to eealieds this by a simplified function of 
the Weibull type, as follows: 


p, . (x) = EXP[ - (x/0 G2) ] 


(:,2) 


In other words, we want oF a. aF > (1,2) to approximate G) Bi, ) 
as closely as possible in a interval 0 = x = Xmax., etuee ial 
is the greatest service period in which we are interested. 

In the first place, we want the relation Tan? ws , 
+ G.) * to be as exact as possible at the very low values of x, say 


0 = x <= 1, because it is most important to have accurate values for 
the probabilities of very short service periods in order that both the 
producer and the consumer may know the true risk of early failure 
for a given rT 

B i . x .B, x ,B, 


=G)' + 





The relation ( 


) 


holds exactly for 








a 2) 
x = 0. If we want it to hold exactly for x = 1, it follows that 
1 1 1 
i a oh ae 
1,2 9, 0 


2 
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From this last equation we are able to solve for the resultant pa- 
rameter Of2), once we know the resultant value of By, ,). Our 
main problem is therefore reduced to one of accurate determination 
of Bi: 2) - 


Determination of the Resultant Weibull Slope *(+,#) 


Before setting up an empirical function for the resultant slope, 
we shall prove the following. 


THEOREM I. The resultant Weibull slope has a value somewhere 
between the two given Weibull slopes, i.e., given two components 
with Weibull slopes B, and B., respectively, the assembly consist- 
ing of these two components has a Weibull slope B (1,2) such that 


PROOF 


Assume 9, < 0, = p®8, (op = 1) 





as closely as possible for all values of x 2 0. 


Putting x = 9, in (1), we obtain 





B 
0, Bia 2) 1 ar 2 
(aaa) 7 +(z) (2) 
0, 
We know that 0;, ,) = 0, =0,. Hence, 0, - =1+€, where e€, 2 0. 
1 2) 
Thus, (2) may be written as 
B B» 
+e) OF 01 oes LALO (3) 
p® pP2 


Furthermore, putting x = 9, in (1) gives 




















Hence, (4) may be written as 
B 
p” (2) (1 + ¢,)° 2) =] + p 1 
Dividing (5) by (3): 


pPt.) = peal + po) 


1 + pB: 


B B 
» 1 1 
or, ath ° Mt = ] 





pP(,2) 1 + pP2 


From (7), it can be seen that 


B, 
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(4) 


(6) 


3 
(A): If B, <B, then x? -! > 1 and, therefore, an — < 1. which 
1+p-2 p (1.2) 


makes B, < B (i 2) - 





B B 3. B 
Furthermore, Sam ES md which makes a) . £ 
2 1 + pB. p® (1 2) pB2 
3 
or > 1 and hence By, .) < B,. Thus, B. < By, 2) < B.,. 
pB( 2) (1,2) ( , ) 1 
1 + pBy pBp 
(B): If By < Bo, then 7. oe <1 and therefore = > 1 which 
+ pBa p® (1,2) 
makes Bi: 2) < B,. 
B B B B 
Furthermore 2— < _7e which makes2®@—*- . P< 4, 
PB. 1+ pBz pB(:,2) pB, 
_ ps, 
or 





pBy, 2) 


<1, and hence B, < By, ,). Thus, B. < By, ,) <B. 
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These deductions indicate that regardless of the relative magnitudes 
of B, and B,, the resultant slope By 2) lies somewhere between 
them. 7 


Q. E. D. 


Since By, 2) falls between B, and B, (Theorem J), it is reason- 
able to express B(, 2) a8 some kind of a weighted average of B, 
and B,. Actual examples of curve fits from past experience have 
indicated that we should express B(i.2) as follows: 


log Ba, 2) = a log B, + (1 - a) log B, (8) 
where Otas i. 
This amounts to the same thing as 
By, .,)=B Bz (9) 


The exponent a is a function of 9, and 0,, and also, in general, a 
function of B, and B,. It should be of such a nature that it will make 
(9) symmetrical in the subscripts 1 and 2. 


If we let + p, a satisfactory formula for a is 
1 





B B, 
on4 (Pe 
re. pB % pB2 
pS + 
For p = 1, this makes (9) to be By, 2) - pt a 4s , the geometric 


mean of B, and B,, which has been Coend to be sioaes by experience 
in the range 0 = x = MAX(Q, ,O,) 


The Fortran Source Program for the IBM 704 

The general synthesis program for Weibull components is based 
on the following equations for the resultant slope and the resultant 
characteristic life of an assembly consisting of two components, (1) 
and (2), with parameters (B, ,0, ) and (B,,0,), respectively: 


RESULTANT SLOPE: By, ,) = B,” B,' * 
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p ” 
+ ; p =>] 
pee a ( 9, 








RESULTANT LIFE: 0.) =—> 





If the assembly in question contains more than two independent com- 
ponents, the above equations are applied repeatedly by treating the 
latest resultant as component (1) and the next component to be con- 
sidered as component (2), until all components are included in the 
final resultant slope and life. 
The life of one of the components should be denoted by unity for 
accurate results. 
The actual Fortran source program looks as follows: 
Cc GENERAL SYNTHESIS PROGRAM FOR WEIBULL COMPONa 
1 FORMAT(120H RESULTANT SLOPE RESUL 
1TANT LIFE NEXT COMPONENT SLOPE AND ITS LIFE 
WRITE OUTPUT TAPE 6,1 
FORMA T(3F24.5) 
FORMA T(4F30.5// ) 
READ INPUT TAPE 7,3,Y,B,X 
Z=1. 
Ww=B 
C=X 
READ INPUT TAPE 7,3,Y,B,X 
WRITE OUTPUT TAPE 6,4,W,C,B,X 
11 R=X/C 
12 BR=(W**((R**B)/(1.+R**B)))*(B**(1./(1.+R**B))) 
12.1 BR1=(W**((R**W)/(1.+R**W)))*(B**(1./(1.+R**W))) 
12.2 BR =SQRTF(BR*BR1) 
XR=(1./((1./(C**W))+(1./(X**B))))**(1./BR) 


ooaouneawt,».* @&© Ww 


= 
oO 





and 
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14 
15 


END 


For illustration we take the following component slopes and lives 
and synthesize them: 


Component Number 


1.00000 
1.66248 
1.34311 
1.15783 
0.99530 
0.99605 


Z=Z+1. 

W=BR 

C=XR 

IF(Z-Y) 9,19,19 
FORMA T(2F30.5//) 
WRITE OUTPUT TAPE 6,18,W,C 


aw h6[}}Se.mlUCU CUD 





AN ACTUAL EXAMPLE 


2.00000 
0.78357 
0.50550 
0.27279 
0.19055 
0.15994 


Weibull Slope 
B, = 


B, = 1.0 


After instructing the IBM 704 to 
source program for the above data we obtain the following output: 


Characteristic Life 


compile and execute the Fortran 


Resultant Slope Resultant Life Next Component Slope And Its Life 





2.00000 
1.00000 
1.00000 
0.50000 
1.00000 


This problem required 1.8 minutes for compiling and 0.1 minute to 
execute, 
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CONC LUSION 


A practical empirical scheme has been presented for quickly 
predicting the resultant Weibull slope and characteristic life of an 
assembly consisting of any number of independent components (in 
series) with known Weibull survivorship functions. This program is 
especially useful in estimating what improvements in the strengths 
of one or more components are required in order to make the relia- 
bility of any final assembly meet specified standards, as well as 
indicating what can actually be expected of some hitherto untested 
assembly made up of tested components. 


APPENDIX 


The Exact Theory 


The exact characteristic life of an assembly consisting of k in- 


dependent components of Weibull slope B,, B,, ..., By, and having 
characteristic lives 0,, 0,, ..., 0, is the positive real root x of the 
equation 
x \By x \Bo x \B 
(o,) + (5) + ee +(3,) - 1 . (10) 


In order to provide a check on the approximate theory presented in 
this report, a Fortran program which solves (10) has been devised. 
This source program is listed below. 


Cc EXACT SYNTHESIS OF WEIBULL COMPONENTS 

1 RATIOF(X,B,A)=(X/A)**B 

2 FORMAT(2F36.5) 

3 FORMAT(2F60.5) 

4 FORMAT(120H COMPON 
1ENT SLOPE COMPONENT LIFE 

5 FORMAT(120H RESULT 
1ANT SLOPE RESULTANT LIFE 


6 X=0. 





a ee ee ee 


g 








A COMPUTER 


7 

8 

9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
33 
34 
35 





PROGRAM 


S=0. 
T=0. 
READ INPUT TAPE 7,2,B,A 
WRITE TAPE 2,B,A 

WRITE OUTPUT TAPE 6,4 
WRITE OUTPUT TAPE 6,3,B,A 
C=.0005 

Y=RATIOF(X,B, A) 

Z=BtY 

S=S+Y 

T=T+Z 

READ INPUT TAPE 7,2,B,A 
WRITE TAPE 2,B,A 

WRITE OUTPUT TAPE 6,3,B,A 
IF(A) 14,22,14 

IF(S-1,) 23,34,34 

X=X+C 

S=0. 

T=0. 

REWIND 2 

READ TAPE 2,B,A 
Y=RATIOF(X,B,A) 

Z=BtYy 

S=S+Y 

T=T+Z 

READ TAPE 2,B,A 

IF(A) 28,22,28 

WRITE OUTPUT TAPE 6,5 
WRITE OUTPUT TAPE 6,3,T,X 
END 
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With the same set of data as in the example on page 39, we ob- 
tain the following printed output from this Fortran program: 


Component Slope Component Life 
1.00000 0.50000 
1.00000 2.00000 
2.00000 1.00009 
1.00000 1.00000 
0.50000 2.00000 
1.00000 1.00000 
0. 0. 


The zeroes at the end of the above table are used to indicate the end 
of the problem. 


Resultant Slope Resultant Life 
0.88753 0.15550 





It can be seen that the previous resultant life 0.15994 as obtained by 
the approximate method is in good agreement with the more exact 
value 0.15550 shown above. The number 0.88753 entitled “Resultant 
Slope” is actually the slope of the tangent at the resultant charac- 
teristic life x = 0.15550 on Weibull paper, and hence, it should not 
be expected to have exactly the same value as the slope of the best 
fitting Weibull line for the assembly. 

The above problem required 1.9 minutes to compile and 2.0 
minutes to execute. 
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Solution of the Orientation Problem 
by Digital Computer Techniques 


By Donald McAllister and Ferdinand Ohnsorg 
Minneapolis-Honeywell Regulator Company 


Abstract 


A solution to the problem of determining the orientation of a 
body in space from the outputs of 3 single axis rotation sensors is 
presented. The solution is accomplished by means of digital tech- 
niques, and suitable algorithms are presented. 


l, INTRODUCTION 


Navigation and guidance of vehicles require knowledge of the 
orientation of the instrument package with respect to inertial space. 
Draper, Wrigley, and Grohe have given a solution to the orientation 
problem by means of gyros mounted on a set of servo driven iso- 
lation gimbals. (1) In some cases, however, it is desirable to de- 
termine vehicle orientation from the outputs of vehicle mounted 
gyros. This article gives a solution to the vehicle orientation prob- 
lem using vehicle mounted gyros and digital computation techniques. 


2. THE ORIENTATION PROBLEM 


Consider a rigid body (Figure 1) rotating with an arbitrary angu- 
lar velocity @ with respect to inertial space. 

Also let 3-single degree of freedom gyros be mounted orthogo- 
nally on the rigid body with input axes along the i, 7, k axes. Let 
the vectors X, 1, y form a non- rotating orthonomal set of unit 
reference vectors in fixed space. The orientation problem then con- 
sists in finding the projections (direction cosines) of the x, BH, = 
vectors along the i, j, k body axes. If the direction cosines of X 
along the 1, j,k axes are Ay» Aas As, respectively the following dif- 
ferential equation between > and w holds, (2) where d is referred to 


the coordinate system of the moving body 
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Y ers = 
Body Fixed 
Space Fixed Axes Gyro Axes 
(non- rotating) (rotating) 
P 
Figure 1. The Orientation Problem 
a — = 
—A=-WxXA 2.1 
dt ( 
with 
w=pi+gqj+rk, 
(2.2) 


The quantities p, q, r are the components of @ along i, j, k respec- 
tively. Similar expressions exist for [i and y. In section 3 we con- 
sider the mechanization of (2.1) by means of an incremental digital 
computer. 


3. INCREMENTAL COMPUTER SOLUTION 
OF THE ORIENTATION PROBLEM 


Let us assume the gyro outputs are sampled at the sampling 
times t,. If 0,, On, Wn are the incremental changes in angle during 
the interval (t,, t,4,), the information shown in equation (3.1) is 
available to us at the time t,,,: 


pdt, on = J qdt, wm = f “rat . a 
t t 
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In this case we restrict ourselves to equally spaced sampling in- 
stants (i.e., sampling occurs at t, = 0, t, = h, --- , t, = nh). 

On the basis of the information in (3.1) we wish to compute, in 
single step fashion, new direction cosines (A,, A,, A, ) at the end of 
the sampling period (as shown in Figure 2). 


| Register | 









































aa) Sr") 
| Syros Sampler > Increment Computer 
Pn, en, Va 
Figure 2, Incremental Direction Cosine Computer 


Before we proceed, let us rewrite equation (2.1) in the more con- 
venient matrix form: 


dy - : o r-q we A, 
-* TA ,: T2#{-F o p > 2&8 A}. (3.2) 
: q-p oO ; As 


We can also define a matrix M of our sampled9,, g,,% from 
(3.1) as 


tn+i oO Vn -Pn 
M,, = f T at -(-4 o On ) ‘ (3.3) 


Ay -9,, oO 


One method of generating suitable relations for the orientation prob- 
lem is to consider the Taylor series of in t. 


= _ h* a* = 
A(tn+,) = . » kl dt* A(t,) , (3.4) 


where h t -t 


n+1 n* 


By repeated use of (3.2) in (3.4) we arrive at 
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Mtns) = [1 + ht +3, (t+ 3) 


(3.5) 
+ h® (T + 2TT+TT + T’) + +++ JA(t ) 
“31 n 
Let us also consider the series 
pm » [rate hte te” 
S Tdt=f] + +3 + 3] becowl, (3.6) 
Or, using (3.3) and (3.6) we arrive at 
h? - 
Mn=hT+ >, T+3, T+--- ; (3.7) 


If both (3.5) and (3.7) are absolutely convergent we can write by 
rearrangement 


Mi Mn 
A(tas,) = [1+ M, + oT + gy A(t) 
(3.8) 


5 
tay (tT - TT)A(t,) + O(h*) . 


Equation (3.8) indicates that the direction cosines , can be 
computed in terms of A(t, ) and Mn, the incremental rotation matrix, 
up to a term of order h* involving the commutator of the matrices T 
NB 
and T 
The mechanization depicted in Figure 2 can thus be carried out 
to order h® by means of the equation (3.8). We will call the term in 
(3.8) of order h® the commutation error 4,at the nth sampling in- 
terval. Hence we have 


— 


h® 
Bn =4g (TT - TTAC,) . (3.9) 


The commutation error 4,, appearing in (3.9) can not be expressed 
in terms of M,, the sampled gyro output, and is inherent in the 
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mechanization of the direction cosine equations (2.1). We will now 
discuss the properties of the commutation error. 


4, COMMUTATION ERROR 


It is well known that rotations do not commute; that is to say 
finite rotations, if treated as vectors, do not obey the vector law of 
addition. Hence, we should not be surprised in this problem that 
digitalization of gyro output gives rise to an intrinsic error A,. 

It is interesting to examine A, to determine the conditions under 
which it is zero. From equation (3.9) it is seen that the commutation 
NB 
error is zero if the commutator T T - T T is zero. 

The quantities T T and T T may be expressed as 


-(rr+qq) qp rp 
tr -( pq -(rr + pp) rq ). 
pr qr - (aq + pp) 
(4.1) 
-(rr+qq) ap pr 
ri =( ap -(rr+pp) rq ) 
rp qr - (4q + pp) 
Hence, commutation error is zero when 
TT-TT=0, (4.2) 


or, equivalently 


Qp-pq=0, 
pr-rp=0, (4.3) 
Qr-rq=0 


The solutions to (4.3) are of the form 


p = p, f(t) 
q, f(t) , (4.4) 
r, f(t) , 


r 


" 


where p,, G,, Tf, are constants and f(t) is an arbitrary function of 
time. 
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The rates characterized by (4.4) are those for which the angular 
velocity of the body @ maintains a fixed direction with respect to 
space. 

We see, then, that commutation error is zero for angular ye- 
locities maintaining fixed orientations with respect to inertial space: 
that is, in the case of planar rotations, finite rotations treated as 
vectors do commute. 

Let us now consider a case where the input rates are given by 


Pp = pp sin at , 


G cos at , (4.5) 


r=0 


From (3,9) and (4.1) the commutation error is 
. ws? 4p-barp-pr A, (ta) 
Re “ja(Pa- ae oO ra- ar) (res) ) - (4.6) 
; ‘. x , 


Or, substituting (4.5) in (4.6) and simplifying 


; h? - A, (ty) 
Bn =e @ Ps oe ( sta) ) ‘ (4,7) 


oO 


Thus the input given by (4.5) gives rise to a computation error whose 
magnitude is essentially a constant per computation cycle. 


5. CONCLUSIONS AND SUMMARY 


a. The orientation problem can be solved by means of digital 
techniques using only the sampled gyro outputs. 


b. An intrinsic error called the commutation error which is of 
the order of the cube of the sampling time interval is present in all 
mechanizations of the direction cosine equations. 


c. The intrinsic commutation error is zero for angular ve- 
locities maintaining fixed directions in inertial space; that is, planar 
rotations. 


d. For sinusoidal angular velocities a non-zero commutation 
error is present. 
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A Method for Obtaining Generalized Transformation 
Operators and An Integral Operator 
By David H. Woodward and Mark D. Epstein * 


Massachusetts Institute of Technology 
Cambridge, Massachusetts 


Introduction 


Operational methods of mathematical analysis have often been 
applied successfully toward the solution of many problems of engi- 
neering and science. In electrical engineering the application of 
operational methods has often been used to simplify the solution of 
differential equations, and in physics the quantum theory has become 
almost synonymous with an operator formalism. Common to all 
possible applications of operational methods, it is clearly necessary 
to possess some knowledge of the appropriate operator. In this 
paper a method for obtaining a generalized type of transformation 
operator will be developed. The application of the method to the 
derivation of a unique integral operator which integrates solely by 
differentiation will be presented. An operator expression for find- 
ing the kth indefinite integral of a function f(x) will be developed. 
Several examples of other applications of the method are also in- 
cluded. 


Generalized Transformation Operator 





Let us first define the differential operator D = =. Before 


treating the general transformation operator, let us first consider 
the example of the common shifting operator, defined as 


[ehD] f(x,) © f(x, + h) (1) 


which is commonly employed inthe calculus of finite differences. (1) 
The operator [e®D] has the property of shifting the value of the 
function f(x) to a new value f(x + h), corresponding to increasing x 
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by a positive unit increment h. This particular operator has been 
employed extensively in the treatment of difference equations, (2) 


Also the operator [e"?} or equivalently, [e®P] where p = = has 


been used in many other operational applications. Oliver Heaviside 
first made use of the property of the operator 


[ehp] f(t) = f(t + h) 


in his operator methods of solution of electrical engineering differ- 
ential equations, (3) 

The operator [e®D] is reallya particular case of a more general 
class of transformation operators to be treated in this paper. To 
help understand the general class of operators, first note that [e"?| 
is defined as 


fehD : = 
=) y Ses n 
ie~l * LZ. =) DM) (2) 


This series form of the operator leads directly to an appreciation 
of the following generalized series form for the transformation 
operator 


[®] 


The constant coefficients A, in (2) are now generalized to be spe- 
cific functions of the variable x. The method to be described next 
is concerned with finding the appropriate coefficients A,(x) corre- 
sponding to a particular operator application. 


[= A,(x)D"] (3) 
n-0 


Method 


The following presentation is intended as a general method for 
deriving transformation operators which can transform any function 
f(x), which is expandable in a Taylor Series, into another desired 
functional form. 

In a fashion similar to expanding a function ina Taylor series, 
the quantity |@|f can be written in an infinite series of powers of de- 
rivatives according to equation (3) in the form 


(O}f = [ 5 An(x)D"|f = A,f(x) + A,f'(x) + Af"(x) + Af" (x) ++: (4) 
r 0 


and each term in (4) can be further expanded in a Maclaurin series 
as 
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2s = , oP a as 
Agf(x) = Ap > oy D"f Ao(f(0) + xf"(0) +>, £"(0) +3) 1”"(0) +--+) 
, S x® oni : ‘ . — 
+ Ajf'(x) = A, Dy DPE = Al (0) + x £"(0) +>, 1"(0) +---) 
, 2 x” ants (5) 
+ A,f"(x) © etn at D*’'*f = Aa( {"(0) + x f”(0) +...) 
o x’ oe 
+ A,f"(x) = A, E75 DPE = As( {"(0) +...) 
. * * * * 


The operator [9] with the appropriate coefficients A,(x), possesses 
the property of being able to transform the function f(x) into another 
desired functional form such as g(x). The problem is to obtain the 
appropriate coefficients A,,(x) corresponding to a transformation to 
g(x). The method by which this may be accomplished procedes in 
the following manner: 


1. Expand the desired functional form g(x) ina Maclaurin series. 


2. Match the g(x) expansion coefficients of derivatives with the 
above operator-function expansion derivatives in equation (5). 


3. From the coefficient values obtained in step 2, write the 
general term for the coefficient An. 


4. Using this expression for A,, write an appropriate form for 


the operator [Qz, ]. 


Two examples of the application of this method will be carried 
through. 


Example 1. 


Consider the example (2) discussed above, where it is de- 
sired to transform f(x) to the form f(x +h), where h is a constant. 
Following the steps outlined above: 


1, Expand g(x), or in this case f(x + h) in Maclaurin series. 


f(x + h) = £(0) + (x + h)f’(0) + {”(0) ["(0) +++. 


(x + h)° (x + h) 
” ae 


2! 





2. Match coefficients of equal derivatives from step 1. with (5). 
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3. Write the general coefficient A,,. 
h®™ 
An - n! 
4. Write a convenient form for the operator. 


is hD 
“4 =[c D"] = 
(Ox+n]=[ 2 2 D®] = [ec] 
Example 2. 


Suppose it is desired to transform a function f(x) into the 
form f(wx), where w is any constant. Then it is a simple matter to 
write out the steps as above: 


1. f(wx) = £(0) + wxf’(0) + a {"(0) + °°: 


2. Match coefficients of derivatives. 





A, = 1 
wk = Aox + A, A, = (w - 1)x 
2 
ius) = men, © ME ORs e° a, = @_) x 
3. The general coefficient is 
{(w - 1)x}" 
An ~ = 
4. Thus 
etx {w= ony | ((w-1)xD 
[@ J [ 2, ene & | = [e 


Many other general transformation operators can be easily derived 
by the above method. The following table lists a few of these oper- 
ators. 





Operator Transformation 
[@, +n] = [e hD] {(x) to f(x + h) 
[O,. J]* erv=?) {(x) to f(wx) 
[Ok ] = [e(x¥-1)xD] f(x) to f(x*) 
[Ok ] = [e(a*-x)D] {(x) to f(a*) 


[O.,bx ] = [e (ax>*-x) D ] f(x) to f(ax>* ) 
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It should be noted at this point that the inverse operators corre- 
sponding to those listed in the table above may also be easily deter- 
mined. Some caution should be exercised in doing this, however, 
because the form of the inverse operator may be different from its 
opposite. For example, while it is true that the inverse of the oper- 
ator [0, + ,] = [e™?] becomes [0,,;,}° =[e®?] corresponding to 
the inverse transformation f(x + h) — f(x), the inverse of the oper- 


: 2% i 
ator (0, ] = [e%-")*P ] becomes [@,,, ]~ =[e'@~*)*] corresponding 


to the inverse transformation f(wx) ~ f(x). 

A brief example of the mechanics of the transformation opera- 
tors will be treated next. Consider the case of the transformation 
f(x) ~ f{(2x), or to be specific, sinh x ~ sinh 2x. 


Example 


Here we apply the second operator listed in the table above, 
with w = 2. Then we have: 


2 3 


{0,, ] sinh x = sinh x + x cosh x + 7H sinh x + i cosh x +... 


2 4 3 5 
, x x x x 
sinh x (1 +9) + Gj ++--) + coshx (x +3) + GF +++) 


= sinh x cosh x + cosh x sinh x 
= 2 sinh x cosh x 


= sinh 2x 


A Unique Integral Operator 





Before discussing the overall significance of the method for ob- 
taining generalized transformation operators, it is perhaps best to 
consider one strange result of this treatment. The result of the ap- 
plication of the method is to obtain an integral operator which finds 
the indefinite integral of a function solely by successive differenti- 
ations. 

The following is a derivation of the integral operator written 
according to the four general steps of the method: 


1, Expand the integral in a Maclaurin series form. 


J-° t(x)dx = x £(0) + 37 1(0) + T f"(0) + z f"(0) + «°- 
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2. Match coefficients of derivatives with equation (5). forn 
this 
A, = x 

x* x 
gy = Aox + A, -» Ay - a (6) 

x* x" x” 
ay = Aogy + Aix + Ao .% A,* 31 Now 
ing 
etc tion 
inde 


3. Write the general coefficient A, . 


; (- 1)®x™* } 
An" Ee Dr 


4. Thus the integral operator is given by ---- 


(- ss" | 


wo J 
-_ > 3 a n © 
l-(Z. Gsm 2! (6 
The 
With the aid of this integral operator, a large class of integrals may alle 
be evaluated in series form without much trouble, Although there for 


are other methods for integrating expandable functions, e.g., inte- 
grating term by term, this operator method has inherently integrated 
the function. In order to integrate a function, it is not necessary to Dis 
expand and then integrate term by term. All that is required is that 
it be possible to differentiate the function. 








To show how the integral operator works, it is probably best to tio 
integrate a simple polynomial. Here however, a slightly more com- has 
plex example of integrating sin x will be given as illustration. Cal 

pay 

Example cal 

the 
[O@;) sin x = x sinx - x cos x - x sin x + x* cos x + x’ sinx - ** sel 
ji 21 31 a1 5! ler 
= Sin x (x - a + r - «++) + cogs x(- a + a ® eset 
= sin’ x + cos x - cOs x 


- cosx + l 


Although most often, the work involved in computing an integral by 
this method is greater, sometimes a quicker, more convenient closed 
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form expression for an integral may be obtained. A fair example of 
this is provided by integrating e*'. 


Example 


fo; e* = e* (x a4 . - Sx’ }=e* 5 (- 1) =a x" 
— © Bae * 8 3. eee 


Now it is possible to extend the above treatment to the case of find- 
ing the kth indefinite integral of a function. By successive applica- 
tions of the method, the integral operator for the 2nd, 3rd, ... , kth 
indefinite integral of a function f(x) may be easily determined. 


eee or, & (1? +1) x” 2, 
Lk f(x, )dx , dx bay in + 2)1 D | f(x) 
. ph, F (-DPm+2)x"” 
{J J, to dx, dx.dx, =[3) + 2 in 3)! D” | f(x) (7) 
x o (- 1)"(n+k-1)x 
LJ, | | f(x, )dx, dx, kl * =W) D" | f(x) ; k2 2 


The use of this operator in finding the kth integral of f(x) will often 
allow faster evaluation and can provide a more instructive closed 
form of expression for the integral. 


Discussion of the Method 

In this paper only a brief summary of the operation and applica- 
tion of the method for finding generalized transformation operators 
has been presented. Except in pointing out the uses of one particular 
case of the generalized transformation operator, the intent of this 
paper has been to present a general treatment of this method appli- 
cable to most transformation situations, Perhaps through the use of 
the type of operators discussed in this paper, it will be possible to 
set up an operator formalism for treatment of transformation prob- 
lems in physical systems. 
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Problems Associated with the Machine Translation 
of Russian Into English 


Principle Investigators: 
Harry H. Josselson 
Arvid W. Jacobson 
Wayne State University 


In translating from Russian into English, or for that matter from 
any language into any other, it is first necessary to discover in each 
language what the various discreet linguistic units are and, second, 
how these units are arranged into larger units in order to carry a 
certain amount of information, After both the source language, the 
language from which we translate, and the target language, the lan- 
guage into which we translate, have been analyzed, it is necessary to 
establish precisely the inter-relationship between both languages on 
all levels—word level, phrase level, sentence level. At this stage it 
becomes necessary to express correspondences between the target 
language and the source language in terms of the former, since it is 
the end output in which we are interested, a finished translation. If 
we reverse the position of the two languages, i.e., make the source 
language the target language, and vice versa, a different set of for- 
mulations will be necessary. If we translate Russian into English, 
we arrive at a series of statements expressing the relationship be- 
tween Russian and English in terms of the final product, which is 
English. In translating from English into Russian, a somewhat dif- 
ferent set of relationships will emerge, if we are interested in ob- 
taining a satisfactory and readable Russian translation. 

Let us illustrate. Here is a simple statement: Russian lacks 
the grammatical category of the article, i.e., Russian does not have 
the three English words THE, A, AN. If we translate from English 
into Russian, we must ascertain under what circumstances the Eng- 
lish article can be omitted in Russian, and if not, what are the con- 
ditions under which the English definite and indefinite articles are 
expressed by differing linguistic structures in Russian. However, 
if we reverse the process and translate from Russian into English, 
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we must ascertain the conditions under which we must add or not aj 
in English, a grammatical category that does not exist in Russian, 
If we are interested in formulating exactly the translatin 
process on the spoken level, the difficulties are, to put it mildly, 
enormous. We must establish the various discreet linguistic units 
on the phonetic (sound alone), phonemic (sound plus meaning), mor. 
phemic and syntactic levels, the inter- relationship of all these levels 
in each language, and then the inter-relationship between each pair 
of languages. We must discover units and establish boundaries, ( 
the written level, we already have well-established units: words, 
phrases, sentences. We have dark masses of type separated by 
white spaces to aid us in establishing units. We have the presence 








of the capitalization device and punctuation mark to aid us in estab. | 


lishing boundaries. Perception of visual symbolism, writing, is 
easier than that of vocal-auditory symbolism—speech. The signals 
of the former are permanently preserved, and can be rescanned if 
necessary, while the latter disappear immediately after they have 
been uttered, unless recorded by some mechanical devices. Both 
types of symbolism, visual, or to be more exact, graphic, and vocal- 
auditory, are of course inter-related and each exhibits its own con- 
plex structure. 

After we have made our separate analyses of both languages and 
established their inter-relationships on the necessary levels from 
the viewpoint of the target language, and after we have done that in 
the symbolism of the natural human language, we have to express 
the whole process in terms of machine language, i.e., in routines 
which a computer can understand. This is the province of the math- 
ematician, logician, and the computer specialist, and the problems 
and complexities involved herein will be dealt with later in this 
paper. Suffice it to say here, however, that the linguist working on 
the preliminary linguistic analyses must be fully aware of the limi- 
tations and restrictions imposed by the computing machine when the 
latter is applied to linguistic analyses. The mathematician, and the 
machine which he uses as a tool, always look for a complete and 
elegant solution: all possible special cases and exceptions must be 
taken into consideration in working out an algorithm, 

We now proceed to a description of our procedures to discover- 
ing the discreet linguistic units of Russian and English and their 
inter-relationship. We are engaged specifically in working out rou- 
tines for translating mathematical texts from Russian into English, 
specifically in the area of partial differential equations. Articles 
from Russian mathematical journals are selected and translations 
of these articles by mathematicians are obtained. We next pro- 
ceeded to edit the English translation with a view of coming as close 
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as possible to the word order in the Russian version. For purposes 
of machine translation, one must also mark separately all those 
items which are ambiguous. For example, in English, TABLE can 
be a noun, a noun-adjunct, as in TABLE LAMP (here TABLE does 
not form a plural, while LAMP does) and a verb: ‘the motion was 
tabled.’ Similarly, in Russian TA3ETN can mean ‘of the newspaper’ 
and ‘newspapers.’ Another example of a more startling ambiguity 
is the following Russian word: HAYAIA, which can mean both ‘she 
began’ and ‘of the principle.’ It is clear, if adequate translation of 
ambiguous items is to be obtained, proper procedures for their 
resolutions must be found. Thus, some words, when treated in 
isolation, can be translated uniquely immediately, while others 
cannot. In addition, a grammatical classification for machine 
translation purposes must of necessity have criteria different from 
that of ordinary analysis. For instance, since a good many pro- 
nouns in Russian are inflected like adjectives, it is better to class 
these items with adjectives. This leaves only the personal, re- 
flexive and KTO - ‘who, which’ and 4TO - ‘what’ in the class of pro- 
nouns, since they are all inflected in a characteristic, easily 
recognizable fashion, differing from all other forms in Russian. 

Thus, each Russian word is marked whether or not it possesses 
characteristics contained in one or more of the following fourteen 
classes: part of speech, declensional class, gender, number, case 
for inflected items or case governed by preceding preposition or 
verb, animation, pronoun type, aspect, reflexivity, tense, person, 
mood, and degree of comparison. Quite obviously, every word 
has a marker in one of the sub-classes of the first class, part of 
Speech, while a noun would be marked, in addition, for declen- 
sional class, number, and case. Animation and gender, in the case 
of nouns, are taken care of in the declensional class system. Thus, 
CTYZEHXT - ‘student’ belongs to a declensional class which is mas- 
culine and animate, in addition to possessing other characteristics, 
such as non-palatalness. Verbs are marked for tense only as 
being either past or non-past. If the verb is imperfective, non- 
past is the present tense, if perfective, non-past is the future. 
Imperfective future is DW, bYMib, etc., plus infinitive. Ambi- 
guities contained within the schemes of declensional paradigms, 
like [TASETH mentioned above, are marked in the declensional 
schemes, while ambiguities arising out of cross-section of form 
Classes, like JAA meaning both ‘they gave’ and ‘of the distance,’ 
will be moved from their original part of speech classes into 
the class of speech labeled ambiguity, after they are discovered. 

By sorting out both Russian and English decks and combining 
them according to the location numbers of the Russian deck, we can 
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construct a glossary, as well as isolate the basic problems involved 
in more-than-word-for-word translation. These are: (1) the prob. 
lem of lexical meaning, unique or multiple, (2) the problem of de. 
letion of certain Russian items in English, (3) the problem of inser- 
tion into English of certain items missing in Russian, and (4) the 
problem of re-arrangement: certain Russian items will have to 
appear in English in a sequence differing from that in Russian in 
order to be understood, 

The glossary must contain in addition to items which are in one. 
to-one correspondence in both languages, with all the problems of 
multiple meaning, solvable by contextual analysis, inherent in 
them, also such occurrences where there is one English word cor- 
responsing to more than one Russian word, and vice versa, like 
TAKMM OBPASOM - ‘thus,’ as well as the problems of multiple mean- 
ings of preposition, idioms and the like. In this connection, the ob- 
vious fact may be pointed out that a glossary compiled on the basis 
of examination of pertinent texts is more up-to-date than one com- 
piled on the basis of examination of existing dictionaries, Russian 
items to be deleted in English, like the emphatic KE, have to be 
suitably marked, We already mentioned above that definite and in- 
definite articles, which are lacking in Russian, will have to be in- 
serted in English wherever absoiutely necessary, although some 
workers in the field would at least like to experiment with the nction 
of having speakers of English read translations of Russian scientific 
texts without any articles contained in them and obtain the readers’ 
reactions. Evolving rules for inserting articles into English turns 
out to be a very complex matter. 

Last, but not least, is the problem of re-arrangement of English 
equivalents of Russian lexical items in a sequence as close to an 
English word order as possible, This involves, of course, pre- 
eminently, syntactic analysis, both within phrases and clauses, as 
well as within sentences. Presumably, after English meanings have 
been assigned to Russian items, they must be further arranged into 
an English sentence. Two approaches are possible: one can analyze 
each Russian sentence, contained between two sentence boundaries, 
by examining each Russian word in turn as it occurs, reaching ten- 
tative conclusions on the basis of this examination and then modify- 
ing these conclusions as further segments of the sentence come into 
view, until a final solution is reached. One could have a series of 
statements, probabilistically expressed or not, on what could or 
could not occur after each item, or what types of patterns are possi- 
ble in Russian sentences, and, what is more important, what the 
English counterparts of these configurations would be. This, in our 
view, as desirable as it is, requires the handling of an enormous 
mass of data, both on the Russian side, as well as on the English 
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side, and in the area of the inter-relationship of both, with a very 
extensive set of logical manipulations. We would, of course, have a 
rather complete picture of the state of affairs, although again not 
quite complete, since certain innovations in sentence structure 
always occur in every language. 

Another approach to analyzing the Russian sentence which may 
be explored, is to assume that each Russian sentence, and English 
for that matter, is an aggregate of groups of words, each of which 
functions as an independent unit in a sentence and can be therefore 
investigated separately. The first problem, then, is to divide the 
sentence into these units, which can be clauses (whether co-ordinate 
or sub-ordinate is immaterial for our purpose) containing in Rus- 
sian, at least, a finite verb or a finite verb-substitute, like predicate 
nouns and adjectives, gerunds, participles, and non-clauses, which 
do not contain a verb or verb substitute. These units must be es- 
tablished both on the basis of boundaries, which are punctuation 
marks or lexical items, like conjunctions, relative pronouns, and 
the like, occurring in the immediate vicinity of these punctuation 
marks, and on the basis of the presence or absence of finite verbs 
or verb substitutes. Each of these units can be analyzed separately 
after it has been isolated. It is hoped that non-clauses, which can be 
parenthetical expressions and the like, can be automatically skipped 
and translated separately, since presumably they will present lesser 
difficulties with respect to re-arrangement when translated into 
English. After the Russian clauses have been analyzed, routines 
will have to be worked out for re-arranging the Russian items when 
they are translated into English, according to a model, which, e.g., 
requires that subject come first, predicate afterwards, that modi- 
fiers always precede the word that is modified, that the items gov- 
erned always follow the governing item, etc. Of course, there is an 
enormous amount of data, first, to be accumulated, next analyzed 
and patterned, and then formulated, It is abundantly clear that, in 
the area of syntax at least, a large number of logical operations will 
have to be performed, It is here that the pay-off in machine trans- 
lation will come and it is here that the mathematician and logician 
can contribute most effectively. 

There are some logical problems involved in unravelling the 
sentence structure. Before these can be adequately considered by 
the logician, a model must be constructed by the linguist. Once this 
has been done, the logician can properly consider the consequences. 
Two models for sentence structure seem tobe in considerable vogue 
among linguists at present, neither of which seems to be suitable for 
machine translation. However, they do serve as a basis for stating 
the logical problem involved. 

a) The Immediate Constituent Model. This is the classic model 
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of sentence structure, familiar to some extent to everyone who has 
diagrammed sentences. To put the model in a different light, we 
may say that between certain words in the sentence a relationship 
exists, as for example, between a noun and its adjective modifier, 
or between a preposition and its object. Two such related words 
may in some sense be regarded as constituting a unit, i.e., an ele. 
ment of some new type in the sentence. The immediate constituent 
model expresses the hypotheses that the structure of the sentence is 
completely defined by the relations existing between adjacent ele- 
ments in the sentence. The sentence structure may, if you wish, be 
represented by a “tree” diagram expressing the successive combi- 
nations of elements into units, as for example, in 


THE GOOD LITTLE BOY HIT THE BAD BIG BOY. 
go ) r \ \ \ \ 


Here, “little” modifies “boy,” “good” modifies “little boy,” and 
the phrase “the good little boy” forms a subject block, while “the 
bad big boy” forms an object block. The phrase “hit the bad big 
boy” is then the predicate, and the whole is a large unit called the 
sentence. 

b) The Transformation Model. This is a relatively new model 
of sentence structure, currently much in vogue among linguists, but 
at present not completely worked out, With the immediate constitu- 
ent model as outlined above, some difficulty is encountered with the 
so-called “discontinuous” constituents as in the sentence 





HE GAVE IT UP. 


Here “gave up” serves as a unit. The transformation model was de- 
signed to handle these, and more serious problems, not satisfacto- 
rily handled by the immediate constituent model, It starts with a set 
of “kernel strings” which are structurally quite simple, and a set of 
“transformations” which operate on the kernel strings to produce 
sentences. The transformations serve to modify word order, insert 
and delete words, change the form of words, and combine several 
structures into one sentence. Now the structure of the sentence is 
specified by listing its set of kernel strings and the set of transfor- 
mations operating on the kernel strings to produce the sentence. 
One of the logical problems involved may be called the decision 
problem. It is the problem of deciding whether a given string of 
symbols is or is not a sentence, The importance of the decision 
procedure is not so much that it enables us to decide if we havea 
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sentence or not, but it gives us implicitly, if it exists, a procedure 
for unravelling the sentence structure. Of course, even with just 
one model, numerous methods for unravelling the sentence structure 
may be available. For example, with the immediate constituent 
model discussed above, it appears to be perfectly feasible to locate 
the boundaries of smaller and smaller units, and thus define the 
structure, by breaking the sentence down rather than building it up. 

As indicated earlier we have selected three articles in Mathe- 
matical text. This was done in order to reduce the problem of mul- 
tiple meaning and perhaps the variability of structural forms. Our 
supposition in the former case has been borne out well. Our analy- 
sis has shown that multiple meaning occurs significantly only in the 
non-inflicted form classes of Russian, particularly the prepositions 
and conjunctions. In the inflected classes the multiple meaning 
problem is confined largely to nouns and adjectives but occurs sur- 
prisingly rarely. 

Another supposition in our approach has thus far been verified: 
general bilingual dictionaries are of little use in automatic machine 
translation. It is necessary to construct one’s own glossary from 
the text itself. In some instances the general dictionary did not even 
include the meaning needed, and in many cases, especially of prepo- 
sitions, the primary meaning in the dictionary was not the primary 
meaning in the text. It appears then that machine translation is 
feasible in relatively narrow subject fields, where special or micro- 
glossaries have been prepared for the given science area and the 
necessary syntactical analysis has been programmed for a computer. 

The syntactic analysis of a sentence proceeds to group words 
into blocks on the basis of agreement or government specified by 
the role of the block in the sentence. The obvious task is to es- 
tablish conditions on such a basis that the computer will be able to 
determine the boundaries of these word blocks. In our approach we 
assume the verb or its substitute to be the pivotal element in the 
sentence. After the principal blocks, including generally nouns or 
nomials, have been identified, we proceed to determine their func- 
tion in the sentence, that is, whether they act, for instance, as the 
subject or the object of the sentence. When the main analysis has 
been completed the final rearrangement, insertions and deletions 
are then performed. 

The machine analysis is envisioned as consisting of several 
passes over the sentence, At the present time we have planned the 
following 10 passes: (1) equation; (2) lookup, missing words, idioms; 
(3) symbols, dashes, numerals; (4) homographs; (5) word blocking; 
(6) inserted structures (subordinate classes); (7) governing modi- 
fiers (participle phrases); (8) gerund packages; (9) relative clauses; 
and (10) main clauses. 
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It may be of interest to mention that we plan to use 108 binary 
bits to describe each Russian word in the machine glossary, Allo | 
the grammatical information, word class, government, agreement | 
codes are included in a fixed bit pattern design. 
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